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A  multilayer  low-energy  electron  diffraction  analysis  of  th&fe!  1 1 1 1  surface  structure  has  determined  four 
interlayer  spacings.  This  analysis  is  a  refinement  of  an  earlier  analysis  that  used  the  same  experimental 
data,  but  considered  only  variations  of  the  first  interlayer  spacing.  The  first  interlayer  spacing  shows  a  small 
increased  contraction,  but  substantial  changes  are  found  in  the  deeper  spacings.  The  optimized  structure  is 
(wit^  the  spacing  between  ith  and  Jrth  atomic  layer)  ^di2  =  (- 16.9+8.0)%,  Arf23=  ( -9  8^3.0)%, 
A ^34  —  (4.2  ±3.6)%,  and  Ad45  =  (  -2.2/±3.6)%,  relative  to  the  bulk  spacing  of  0.827  A. 


I 


■  ■'<4 


I.  INTRODUCTION 

The  occurrence  of  multilayer  relaxation  of  the  outer 
layers  of  metal  crystals  is  now  well  established.1  About  ten 
cases  have  been  studied  by  low-energy  electron  diffraction 
(LEED)  and  ion  scattering.  In  particular,  the  study  of  six 
surfaces  of  Fe  has  shown  that  the  relative  contraction  of  the 
first  layer  spacing  d\2  increases  smoothly  with  the  openness 
or  roughness  of  the  surface  and  that  relaxations  go  several 
layers  deep.2  The  original  work  in  the  Fell  11)  surface,3 
which  was  studied  along  with  the  other  low-index  surfaces 
(100)  and  {110}  and  varied  only  dl2,  found  the  surprisingly 
large  contraction  of  ( 1 5.4  ±3 )%.  However,  now  that  the 
dependence  on  openness  is  known,  this  result  for  the  (111) 
surface,  a  rather  open  surface,  no  longer  appears  as 
anomalous.  What  appears  to  be  missing  is  a  study  of  the 
changes  in  deeper  interlayer  spacings,  particularly  because 
open  (i.e.,  loosely  packed)  metal  surfaces  are  now  known  to 
exhibit  large  multilayer  relaxations.2  Of  interest  is  also  the 
question  about  what  change  would  be  found  in  the  magni¬ 
tude  of  d\2  when  multilayer  relaxation  would  be  allowed. 
Thus,  we  have  made  a  new  study  of  the  Fe  { 1 1 1 J  surface 
taking  into  account  changes  in  the  first  four  interlayer  spac¬ 
ings.  The  results  show  -that  such  an  analysis  is  necessary, 
since  substantial  relaxation  of  deeper  layers  occur.  Also, 
since  the  early  work  was  done,  the  theory  of  relaxation  of 
metal  surfaces  has  advanced  to  the  status  of  first-principles 
calculations4  which  predict  the  sequence  of  multilayer  relax¬ 
ations.  These  predictions  are  based  on  total-energy  self- 
consistent  calculations  which  relax  layer  positions  as  well  as 
the  electron  distribution.  Hence  the  present  relaxation 
results  provide  a  valuable  test  of  those  theories.5 


II.  ANALYSIS 

The  experimental  data  base  of  14  nondegenerate  LEED 
spectra  at  two  angles  of  incidence  used  in  the  previous 
Fell  11)  analysis3  has  been  reanalyzed  using  as  adjustable 
parameters  the  first  four  interlayer  spacings,  d\2,  d2 3,  dyt, 
and  d* 5  and  the  real  part  of  the  inner  potential  F0.  Theoret¬ 


ical  spectra  were  calculated  with  the  same  computer  pro-  / 
grams  and  scattering  potential  used  for  previous  iron  stud¬ 
ies.6  The  imaginary  part  of  the  inner  potential  was  fixed  at 
4  eV,  and  up  to  55  beams  were  used  to  represent  the  wave 
function  between  layers.  The  agreement  between  experi¬ 
mental  and  calculated  spectra  was  evaluated  quantitatively 
with  the  numerical  r  factor  of  Zanazzi  and  Jona.7 

The  method  used  to  optimize  the  structural  and  nonstruc- 
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FIG.  1L  Experimental  and  theoretical  LEED  spectra  for  Fell  11), 
for  10,  TO,  and  11  beams  at  9-0°.  EXPT.  -  experiment;  F.R. 
-fully  relaxed  structure;  T.L. -top-layer  relaxation  only. 
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tural  variables  was  described  previously.8  Briefly,  the  op¬ 
timization  was  done  as  follows.  (1)  A  series  of  calculations 
involving  independent  variations  of  all  the  structural  param¬ 
eters  by  relatively  large  amounts  were  done  in  order  to  get  a 
rough  approximation  of  the  optimum  structure.  For  all  cal¬ 
culations,  the  value  of  the  real  part  of  the  inner  potential  VQ 
was  allowed  to  vary  independently  and  a  minimum  r  factor 
was  found.  (2)  All  but  two  of  the  structural  parameters 
were  fixed  at  the  cAide  “best”  values  found  in  (1)  and  the 
remaining  two  structural  variables  were  allowed  to  vary  in¬ 
dependently  until  a  minimum  of  the  reliability  factor  was 
found.  With  the  refined  values  of  these  two  parameters 
thus  obtained,  a  new  pair  of  variables  was  chosen  and  the 
process  repeated.  The  pairs  considered  were 
(<*23. <*34),  and  (di4,d45)\  when  the  consistency  among  op¬ 
timum  values  determined  for  members  of  different  pairs 
was  within  the  estimated  experimental  uncertainties 
(3%-4%  of  an  interlayer  spacing),  the  analysis  was  ter¬ 
minated.  In  addition,  variation  of  the  angle  of  incidence 
away  from  the  nominal  experimental  value  of  9°  (measured 
on  the  sample  goniometer)  was  considered  and  an  optimum 
value  of  8°  was  obtained.  The  results  of  the  analysis  were 
as  follows: 


FIG.  2.  Experimental  and  theoretical  LEED  spectra  for  Fell  11), 
for  20  and  22  beams  at  0  —  0°.  EXPT.  —  experiment;  F.R.- fully  re¬ 
laxed  structure;  T.L.  —  top-layer  relaxation  only. 


<*12  =  0.69  ±0.025  A  [( —  16.9  ±3.0)%  contraction  relative  to  bulk  value  of  0.827  A]  , 
<*23  =  0.75  ±0.025  A  [(  —  9.8  ±3.0)%  contraction]  , 

<*24  =  0.86  ±0.03  A  [( +4.2  ±3.6)%  expansion]  , 

<*45  =  0.81  ±0.03  A  [(-2.2  ±3.6)%  contraction]  , 

Vo,-  -  12.6  ±  1.0  eV  (real  part  of  the  inner  potential)  , 
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FIG.  3.  Experimental  and  theoretical  LEED  spectra  for  Fell  11), 
for  00,  10,  and  ll  beams  at  0-8°,  <3-0°.  EXPT.  -  experiment; 
F.R.  -fully  relaxed  structure;  T.L. -top-layer  relaxation  only. 
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FIG.  4.  Experimental  and  theoretical  LEED  spectra  for  Fell  11), 
for  12,  21,  and  10  beams  at  0-8°,  <3-0°.  EXPT.  -  experiment; 
F.R. -fully  relaxed  structure;  T.L. -top-layer  relaxation  only. 
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FIG.  5.  Experimental  and  theoretical  LEED  spectra  for  Fell  11), 
for  31,  20,  and  21  beams  at  9  —  8°,  EXPT.  -  experiment; 

F.R.  -  fully  relaxed  structure;  T.L.  -  top-layer  relaxation  only. 


(/■min  is  the  minimum  of  the  Zanazzi-Jona  r  factor  for  the 
data  base  of  14  beams;  the  value  obtained  when  only  varia¬ 
tion  of  the  first  interlayer  spacing  is  considered  is  0.156). 

Figures  1-5  show  the  comparison  between  experimental 
and  theoretical  LEED  spectra  for  both  the  fully  relaxed 
model  and  the  optimum  model  obtained  when  only  top- 
layer  relaxation  was  considered.3  Figure  6  shows  a  plot  of 
top-layer  surface  relaxation  versus  surface  roughness  or 
packing  fraction  (top  scale)  for  six  iron  surfaces.2  In  Fig.  6, 
the  cross  and  the  circle  for  the  (111)  surface  indicate  the  op¬ 
timum  values  for  top-layer  relaxation  found  when  consider¬ 
ing  top-layer  relaxation  only  and  multilayer  relaxation, 
respectively.  The  squares  give  the  relative  total  relaxation9 
of  the  first  layer  for  the  surfaces  with  multilayer  relaxation. 


m.  CONCLUSION 

The  multilayer  analysis  has  increased  the  relative  contrac¬ 
tion  of  d\2  to  16.9%  and  revealed  a  large  contraction  of 
which,  combined  with  changes  in  d 34  and  d45,  gives  a  total 
relaxation  Ad  of  the  first  layer  with  respect  to  bulk  of  -0.20 
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FIG.  6.  Top-layer  relaxation  of  Fe  surfaces  relative  to  bulk  spac¬ 
ing  in  percent  vs  surface  roughness  (inverse  of  packing  fraction). 
The  cross  is  the  result  obtained  when  the  relaxation  of  the  top  layer 
is  the  only  relaxation  allowed  (Ref.  3);  the  circle  is  the  result  ob¬ 
tained  with  the  refinement  done  in  the  present  work.  The  squares 
give  Ad/d12,  the  total  relaxation  of  the  top  layer  (Ref.  9)  relative  to 
bulk  spacing  in  percent. 
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A,  the  largest  valu^  of  the  six  iror^  surfaces  investigated  to 
date  (Ad-  -0.08  A  (211),  -0.07  A  (310),  -0.15  A  (210)), 
and  about  the  same  relative  total  contraction  as  (210) 
(Ad/d12--6.8%  (211),  -7.5%  (310),  -24%  (210),  -24% 
(111)).  The  sequence  of  changes  in  layer  spacings  (contrac¬ 
tion,  contraction,  expansion,  contraction)  is  the  same  as  for 
the  (210)  surface,  whereas  the  (211)  and  (310)  surfaces 
show  strict  alternation  of  sign.  Determination  of  the  signs 
and  magnitudes  of  relaxations  on  these  various  surfaces 
provides  good  quantitative  data  to  test  metal  surface  theory. 
The  new  multilayer  parameters  for  Fell  11)  offer  the  unique 
combination  of  a  surface  that  is  both  strongly  relaxed  and 
highly  symmetric,  without  the  complication  of  parallel  relax¬ 
ation.10 
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